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Abstract 

In the analysis of thermal infrared data of asteroids by means of thermophysical models (TPMs) it is a common prac¬ 
tice to neglect the uncertainty of the shape model and the rotational state, which are taken as an input for the model. 
Here, we present a novel method of investigating the importance of the shape model and the pole orientation uncertain¬ 
ties in the thermophysical modeling - the varied shape TPM (VS-TPM). Our method uses optical photometric data to 
generate various shape models that map the uncertainty in the shape and the rotational state. The TPM procedure is 
then run for all these shape models. We apply the implementation of the classical TPM as well as our VS-TPM to the 
convex shape models of several asteroids together with their thermal infrared data acquired by the NASA’s Wide-held 
Infrared Survey Explorer (WISE) and compare the results. These show that the uncertainties of the shape model and 
the pole orientation can be very important (e.g., for the determination of the thermal inertia) and should be considered 
in the thermophysical analyses. We present thermophysical properties for six asteroids - (624) Hektor, (771) Libera, 
(1036) Ganymed, (1472) Muonio, (1627) Ivar, and (2606) Odessa. 
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1. Introduction 


The physical characterization of asteroids has seen an en ormous boost in recent years thanks to th e data of the 


NASA Wide-held Infrared Survey Explorer (WISE) mission (Wrieht et al. 201C; Mainzer et al. 

2011a 

). 

In particu- 

lar, sizes and albedos have been determined from simple thermal models (e.g., NEATM of 

than 150,000 asteroids and for different populations thereof, including main belt asteroids (' 

Harris 

1998 

) for more 

VlBAs) - 

Masiero et al. 


(2011 


2012), Hildas- 


Gray et al 


(12011 


Gray et al 


(l2012al) . near-Earth objects (NEOs) - 


Mainzer et al. 


(12014 


201 Ibl) . and Trojans - 


2012d) . This resulted in a database of impressive quality in terms of number of observed bodies 
and sensitivity as compared to previous surveys: ~2200 asteroid albedos and sizes from the SIMPS based on IRAS 

- II I -^-^1- 

observations (ITedesco et al.ll2002^ . or ~5 000 albedos and sizes based on AKARI data (lUsui et al.ll201 II) . 
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Simple thermal models assume spherical, non-rotating asteroids. When shape models, spin vectors, and other 
physical parameters are known, it is possible to apply more sophisticated methods to infer thermophysical properties 
such as thermal inertia and/or surface roughness. So far, such thermophysical parameters were derived for less than 
about 50 asteroids due to the lack of both physical properties (i.e., shape models, spin axis orientations) as well as 
available thermal infrared observations. However, the number of shape models has grown dramatically in the past 
decade and the newly obtained thermal infrared data from WISE now open the opportunity to greatly increase the 
number of thermal inertia determinations. 

Thermal inertia, defined by E = (pKC)^^^, where p is the density of the surface regolith, k its thermal conductivity, 
and C its heat capacity, measures the resistance of a material to temperature change, and thus controls the temperature 
distribution of the surface of an atmosphere-less body. As it affects the symmetry of the temperature distribution on 
asteroids, the thermal inertia controls the strength of the Yarkovsky effect, which is the rate of change in the s emi- 


major axis of the orbit of an asteroid (da/dt) due to the recoil force of the thermal photons jBottke et al 


20061 for 


instance). Thermal inertia is also a sensitive indicator of the nature of the surface regolith as its value is affected 


e.g., 


Mellon et al. 

2000 

Jakoskvl 

1986) 


for future landing and/or sample-ret urn missions (such as OSIRIS-REx and 


missions to carbonaceous a steroids. 


Lauretta et al. 

2 f )12 

Okada et al. 

2014) 


values of both k and T. See 




Vernazza et al 


Hayabusa-2, two different sample return 


20141) . Thermal inertia is strongly affected by 


the porosity of the material dZimbelman|2986^^or a given surface composition, the higher the porosity, the lower the 


( 12012 h for a discussion of the effect of porosity on asteroids surfaces. 


In order to derive the thermal inertia and oth er physical parameters of as teroids, such as the diameter and the 


albedo, a thermophysical model (hereafter TPM, 


Lagerros 


1996 


1997 


19981) is typically used to analyze thermal 


infrared data. A TPM calculates thermal infrared fluxes given a set of physical parameters (size D, thermal inertia 
E, Bond albedo A, surface roughness 0) whose values are adjusted to provide the best fit between the model and the 
observed fluxes, by minimizing a figure of merit (chi-square). Classically, a TPM is used with an a-priori knowledge 
of the shape and the rotational state of the asteroid, which are take n as fixed quantities. Typically, shapes are based 


on radar imaging (e.g., asteroids 2010 EV 5 , or (101955) Bennu, 


All-Lagoa et al 


2014 


Emery et al 


convex inversion of photom etric ligthcurves (e.g., asteroids (25143) Itokawa, or (1620) Geographos, 


2014) or on 


Muller et al 


2014 


Rozitis & Greenll2014t). 


So far, shape uncertainties have never been properly accounted for in TPM analyses, but there is growing evidence 
that large chi-squared values are obtained when the s hape model used f or the thermophysical modeling differs sig¬ 


nificantly from the asteroids true shape. In particular. 


Rozitis & GreenI (120141) have shown that if the length of the 


shape models rotation axis, as determined from radar observation, is not properly estimated, this can lead to large chi- 


2014 


Muller etal .1120121). 


squared values and bias t he value of the thermal inertia (see also the case for the NEA (101955) Bennu bv lEmerv et al 

|[ 


Here we study how the uncertainties of the shape model and the pole orientation influence the TPM results; we 
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perform TPM analysis of WISE thermal infrared data for a selected number of asteroids t hat have known shap es and 

itQ 


spin poles from the Database of Asteroid Models from Inversion Techniques (DAMIl 


Durech et al 


2010) . First, 


using new optical lightcurve data, we derive revised shape models that we use for TPM analysis. We show that the 
goodness of fit of the TPM (e.g.,;^f^) changes when the revised shape is used instead of the DAMIT one, indicating 
that WISE data are sensitive to shape features. We also note that the revised shape models can, in general, improve 
the TPM solution. 

Next, we introduce a novel method, that we call the varied shape TPM (or VS-TPM for short), that allows us to 
analyze the stability of the TPM solution against the variations of the shape model and the pole orientation (Sect. |4]). 
We show that our method also allows us to find a shape model that improve significantly the fit of thermal infrared 
data compared to a classical TPM based on fixed shape and pole orientation derived uniquely from optical light curve 
inversion (nominal shape). 

The VS-TPM c onsists of the following s teps: w e bootstrap the opt i cal ph otometric data and use the technique of 


convex inversion of 


Kaasalainen & Tomnal (1200 ih : 


Kaasalainen et al 


(120011) to determine a set of slightly different 


shape solutions that fit the available optical disk-integrated photometry equally well. For each shape and pole solution 
(varied shapes), we use a TPM to analyze thermal infrared data (e.g. WISE, IRAS) and we notice that the goodness 
of the fit to thermal infrared data can be significantly shape model dependent, and in most of the cases the VS-TPM 
allows one to find a solution that fits the thermal infrared data better than the original shape solution. 

In this paper, we first describe and apply in Sect.[3]the classical TPM scheme with an a priori convex shape model 
to the WISE data for nine selected asteroids. In Sect.|4] we present the VS-TPM method and use it to show how the 
uncertainties in the shape model and its pole orientation influence the thermophysical fit. We discuss and conclude 
our work in Sects. |5]and|6] Application of the varied shape TPM scheme to a few hundred of main-belt asteroids will 
then be a subject of our forthcoming work. 


2. Data 

2.1. Visible light, disk-integrated photometry 


*littp: //astro. troja .m££. cuni. cz/projects/asteroids3D 


3 














Table 1: List of asteroid models derived by the lightcurve inversion method that we use for the TPM modeling. For each asteroid, the table gives 
the ecliptic coordinates A and /3 of the pole solutions, the sidereal rotational period P, the number of dense lightcurves A^ic observed during A^app 
apparitions, the number of spai'se data points from USNO-Flagstalf A ^689 Catalina Sky Survey A^ 703 , and the reference. 


Asteroid 

M 

[deg] 

ySi 

[deg] 

^2 

[deg] 

(^2 

[deg] 

P 

[hours] 

N\, 

^app 

Nm 

Njoi 

Reference 

624 Hektor 

331 

-32 



6.92051 

17 

8 



Kaasalainen et al. (2QQ2a) 


333 

-32 



6.920509 

19 

9 

201 

56 

This work 

771 Libera 

64 

-78 



5.89042 

20 

5 



Marciniak et al. ('20091 

832 Karin 

242 

46 

59 

44 

18.3512 

13 

6 

84 

39 

Hanus et al. ("20111 

1036 Ganymed 

214 

-73 



10.313 

21 

1 



Kaasalainen et al. 12002bl 


190 

-78 



10.31284 

177 

4 

155 

20 

This work 

1472 Muonio 

42 

62 

249 

61 

8.70543 

6 

1 

99 

93 

Hanus et al. 1201361 

1627 Ivar 

338 

40 



4.79517 

56 

4 



Kaasalainen et al. 120041 


334 

39 



4.79517 

83 

4 

68 

152 

This work 

1865 Cerberus 

292 

-72 



6.803284 

28 

8 



Durech et al. 120121 


311 

-78 



6.803286 

47 

8 


62 

This work 

1980 Tezcatlipoca 

334 

-65 



7.25226 

48 

5 



Kaasalainen et al. 120041 


324 

-69 



7.25226 

49 

6 

29 

35 

This work 

2606 Odessa 

25 

-81 

283 

-88 

8.2444 

3 

1 

25 

129 

Hanus et al. 1201361 
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Table 2; Optical data used for the shape model determinations. 


Asteroid Date N\c 


624 Hektor 

1957 4- 1957 5 

4 


1965-02-04.3 

1 


1967-03-07.4 

1 


1968 4- 1968 5 

2 


1977 2- 1977 2 

2 


1984 10- 1984 10 

3 


1990 3 - 1990 3 

2 


19914- 19914 

2 


2008 10-2008 10 

2 

771 Libera 

1984 5- 1984 5 

3 


1999 9 - 1999 9 

2 


1999 9 - 1999 9 

2 


2005 2 - 2005 3 

3 


2006 5 - 2006 6 

4 


2008 10-2009 3 

6 

832 Karin 

1984 10- 1984 10 

2 


2003 8 - 2003 9 

8 


2004 9 - 2004 9 

3 

1036 Ganymed 

1985 7- 1985 11 

6 


1985 7- 1985 12 

25 


1989 4- 1989 7 

11 


2008 12-2009 4 

24 


2011 5-2012 1 

8 


2011 5-2011 12 

103 

1472 Muonio 

2008 9 - 2008 9 

3 


2008 10-2008 10 

3 

1627 Ivar 

1985 6- 1985 10 

26 


1990 5- 1990 8 

18 


1990 5- 1990 5 

2 


1990 5- 1990 5 

2 


1995 2 - 1995 3 

6 


Observer 


Dunlu£_&_Gehr^ (1969) 
^unlag,&^ehrels (]%9) 
_Dunla£&_Gehrels (]%9) 
Dunlu£_&_Gehr^ (]%9) 
Ifetmann_&Cruikshank (1^7^ 
^etalet^ (]^9A) 
^ahl^ren_et^ (J39V) 
ffainaut;^ouelle_et^ (19951 
^teghens(2009^ 

Binzel (19871 
Marciniak et al. (20091 
Warner (20001 
Marciniaket^ (2009) 
Marciniak et al. (20091 
Marciniak et al. (20091 
Binzel (19871 
Yoshidaet^ (200A) 

Ito & Yoshida (20071 
Lupishko et al. (19871 
Hahn et al. (19891 
Chemova_e^^ (19951 
Skiff et al. (20121 
Pilcher_e^^ (201^ 
Velichko et al. (20131 
Stephens (2009bl 
Hanus et al. (2013bl 
Hahn etal. (19891 
Chemovae^^ (19951 
Velichko_et^ (1^90) 
Iloffinann_&^e^er (1V90) 
Pravec et al. (19961 


Observatory (MFC code) 


Lowell Observatory 
Multiple observatories 
Multiple observatories 


Ondfejov Observatory (557) 
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Table 2; continued. 


Asteroid Date N\c Observer Observatory (MFC code) 


2008 9 - 2009 2 29 

1865 Cerberus 1980 11 - 1980 11 2 

1989 11 - 1989 11 2 

1998 10- 1998 10 2 

1999- 09-25.0 1 

1999 11 - 1999 11 2 

2000- 7-09.1 1 

2008 9 - 2008 11 25 

2008 10-2008 11 5 

2009 9-2009 10 4 

2010-08-06.9 1 

1980 Tezcatlipoca 1988 6 - 1988 6 2 

1992-05-22.0 1 

1995-10-27.1 1 

1996 2 - 1996 2 3 

1997 6- 1997 1 41 

2009-08-20.1 1 

2606 Odessa 2008 2-2008 2 2 

2008-03-1.5 1 


Skiff etal. (20121 
Harris & Young (1989') 
jMsniewski_et^ ( 1997 ) 
^mieczk^et^ (199^ 
^zab6_et^ (20^ 
^Dutech_et^ (201^ 
^zab6,et^ (2001) 
Durech_et^ ( 2012 ) 
Skiff et al. (20121 
Durech,et^ (201^ 
Durech,et_^ (201^ 
JMsniewski_et^ ( 1097 ) 
ICaasalainen_e^^ (200£) 
ICaasalainen,e^^ (200^ 
ICaasalainen,e^^ (20QA) 
Kaasalainen ef al. (20041 


Skiff etal. 120121 
Hig2ins.et^ (200^ 
JHanus_e^^ (2013^ 


Lowell Observatory 


We make use of optical photometric data for two reasons; (i) to revise/improve several shape models, and (ii) to 
bootstrap the photometry to derive various shape models for the VS-TPM. 

We use two different types of reflected disk-integrated photometry; (i) dense-in-time photometry, which is typ¬ 
ically acquired by individual observers and densely covers a time interval of several hours, and (ii) sparse-in-time 
photometry, which is a usual by-product of astrometric surveys and consists of a few hundred individual calibrated 
measurements during ~15 years. 

Shape models adopted from the literature (stored in the DAMIT) are usually based on the dense-in-time pho- 
tometry. The dense -in-time photometry is from two main sources; (i) the Asteroid Photometric Catalogue (APC 


Piironen et al. 


20011), and (ii) the data from individual observers provided by the Minor Planet Centeij^ in the Aster- 




^http://asteroid.astro.helsinki. £i/ 

'^http : //WWW .minorplanetcenter. net/light_curve2/light_curve. php 
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oid Lightcurve Data Exchange Format ('ALCDEIri IWarner et al.ll2009h . Several observers send us data directly or 
on request. Our revised shape model determinations are based on the combined dense- and sparse-in-time data sets. 
We download the sparse-in-time photometric data (typically a cquired by a s trome tric surveys) from the AstDyS site 


(Asteroids - Dynamic 


ic Sitjf ) 


and process them similarly as in 


Hanus et al 


(1201 ih : we compute the geometry of the 


observations, light-time correct the epochs, transform magnitudes to intensities, and exclude clear outliers. We use the 
sparse photome tric data from the U SNO-Flagstaff station (lAU code 689) and the Catalina Sky Survey Observatory 
(lAU code 703, Larson et alJl2003 ). 


In Tab. [T] we list for each studied asteroid information about the visible photometry used for the shape model 
determination, namely the number of dense-in-time lightcurves, the number of apparitions covered by dense-in-time 
observations and the number of sparse-in-time measurements from both astrometric surveys. The references to the 
photometric observations are presented in Tab. 

2.2. WISE thermal infrared fluxes, disk-integrated photometry 

We make use of the thermal infrared data of asteroids acquired by the WISE satellite, in part icular the results of the 
NEOWISE project, which focuses on the solar system bodies (see, e.g.. 


Mainzer et al.ll201 lal) . The thermal infrared 

data are downloaded from the WISE All-Sky Single Exposure Lib Working Database via the IRSA/IPAC archive^. 

In this study, we consider only thermal infrared data from filters W3 and W4 (isophotal wavelengths at 12 and 
22 pm) from the fully cryogenic phase of the mission. While W3 and W4 data are thermal-emission dominated, the 
fluxes in filters W1 and W2 (isophotal wavelength at 3.4 and 4.6 pm) usually at least partially consist of reflected 
sunlight, which cannot be properly modeled in our purely thermophysical model and these filters are therefore not 
considered here. 

The data selection and suitability criteria applied in this work follow those of 


(34184 3) 2 008 EVs. In turn, t hese are based on a combination of criteria from 


2011 

), and 

Grav et al. 

(2012a) 


AlELagoaet^ 


Mainzer et al.l (1201 id) . 


(2014) for asteroid 


Masiero et al 


(l2012al) : we implement the correction to the red and blue calibrator discrepancy in W3 and W4 


(Cutri et al. 2012), and we use a cone search radius of 1" centered on the MFC ephemeris of the object in our queries. 
We only consider data with artefact flags p, P, and 0, and quality flags A, B, and C. A “0” artifact-flag e ntry indicates 


no arti fact detection, whereas p and P indicate possible contamination by a latent image. Nonetheless, 


Masiero et al 


(1201 ll) found the pipeline criteria for P and p flags to be overly conservative, in particular, that one safely retrieves 20% 


more data by allowing P-flagged. The quality flags A, B, and C, correspond to the following signal-to-noise ratios 
(S/N): S/N > 10 (A), 3 < S/N < 10 (B), and 2 < S/N < 3 (C). We require the IRSA/IPAC modified Julian date 
to be within four seconds of the time specified by the MFC. A positive match from the WISE Source Catalog within 
6" of any MPC-reported detection indicates that there is an inertial source at a distance smaller than the point-spread 


'^http : //WWW . minorplanet. in£o/alcde£. html 
^http://hamilton.dm.unipi.it/ 

%ttp://irsa.ipac.caltech.edu/Missions/wise.html 
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function width of band W1. We consider that these data are contaminated if the inertial source fluxes are greater than 
5% of the asteroid flux and we remove them. Finally, we do not use data that are partially saturated to any extent. 


Additionally, we And indications that the error bars of the WISE data may be underestimated, namely the fact 
that those measurements taken 11 seconds apart are not usually compatible within their error bars. These “double" 
measurements occur when the asteroid appears in the overlapping area of two consecutive frames. We identify ~400 
such double detections in the WISE data set for asteroids with convex shape models, and we And that the reported 
uncertainties of the thermal fluxes are underestimated by a factor of ~1.4 for W3 data and of ~1.3 for W4 data 
(described in more details in Appendix A). The change in the thermal flux of an asteroid during this short time 
interval is much smaller than the relative error of the fluxes and it cannot be caused by the orbital and rotational 
evolution either, therefore we decide to increase the error bars of the data by these factors. 

Eurthemore, 


Jarrett et al. 


(1201 Ih presented a study of the WISE-Spitzer flux cross-calibration of a number of 


calibration stars and one galaxy situated near the poles of the ecliptic and found that (i) the photometry for individual 
objects is stable for the whole cryogenic phase within less than 1%, but (ii) there is an rms scatter around the zero level 
of 4.5% and 5.7% of the WISE zero magnitudes in filters W3 and W4 when several objects are examined (the offset is 
different for each individual star). Accounting for these uncertainties in the case of asteroids is questionable, because 
the cross-calibration is studied for stars. Asteroids, however, have important differences compared to said calibration 
stars: they have widely different surface temperatures ~4000 vs. ~300 K; asteroids can present rotational variability 
of a factor of two or more in their fluxes due to their shapes; also the typ i cal ma gnitudes of asteroids in Alters W3 


and W4 are much brighter than those of the stars examined by 


Jarrett et al. 


(120111) . In this context, it is reasonable to 


expect an unknown and different offset for each asteroid. It may even change within the dataset of the same object, 
due to its intrinsic time variability. Thus it can only be accounted for in our modeling by a random relative offset of 
the fluxes between the W3 and W4 bands. Our method to deal with this source of error and its effect on the physical 
parameters of asteroids derived my means of TPMs is described in Appendix A. In short, we And that the error of the 
cross-calibration uncertainty does not affect the results significantly. 


3. Thermophysical modeling with fixed shape models 

In order to study the importance of shape uncertainties in the thermophysical modeling, first we study nine selected 
asteroids for various reasons, for example, high quality shapes or wise data, or recent additional optical data are 
available. We then show that by applying the TPM on a revision of the shape model derived from the inversion of 
new lightcurve data in addition to the one present in the DAMIT, we obtain a different goodness-of-fit. In general, our 
results improved, but the revised best-fitting parameter values are consistent within the error bars. On the other hand, 
still stay above 1 in most cases. 
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Figure i: Dependence of the values of the TPM fits on the themial inertia F for a rotational phase = 0° and five different surface roughness 
values 6 for the revised shape model of asteroid (1627) Ivar. 
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3.1. TPM method 


We use a thermophysical model imp l ementation of 


developed bv iLagerrosI (11996 


1997 


19981); 


Spencer et al 


19891) : 


Spenceij(ll990l) : 


Delbo’ et al. ( 2007 ): Delbo’ (2004) that is based on TPM 


Emery et al 


(Il998l) . A TPM allows 


thermal infrared fluxes to be calculated at different wavelengths and at a number of epochs taking into account the 
shape of an asteroid, its spatial orientation, and a number of physical param eters such as the size of the body, the 


albedo A, the macroscopic surface roughness 6 (Hapke’s mean surface slope. 


Hapke 


19841) . and the thermal inertia 


r. The values of the parameters are determined by minimizing the difference between the observed fluxes and 
the modeled fluxes s^Fi, where we consider the scale factor s for the asteroid size, and i corresponds to individual 
observations. To find the optimal set of parameter values, we minimize the metric 






(s^Fi-0 


( 1 ) 


where cr, represent the errors of fluxes /). 

The shape is represented by a convex polyhedron with triangular facets (models f rom the DAMIT database) o r 
by a set of surfaces and normals (for revised models, the so called Gaussian image, see 


Kaasalainen & Torpp^l^OOll) 


According to our tests, both representations produce similar fluxes and thus are equivalent. In the case of new models, 
we favor the latter representation, which is a natural outcome of the convex inversion, because it is less computationally 
demanding. The shape has an initial size, and the parameter s, adjusted, to FIT the data, is a factor that either scales 
linearly all vectors of the vertices of the polyhedron, or quadratically all the surfaces of the Gaussian image. Given 
the best-fitting value of s, we determine the volume equivalent diameter D (i.e., the diameter of a sphere with the 
same volume as the scaled shape model), which is a quantity we present throughout this work and that we use to 
determine the body’s albedo. For the WISE bands, the bolometric emissivity can be set to a constant value of 0.9, 
that is the typical average spectral emissivity in the range of thermal infrared wavelengths used in thermal models. 
At the relevant temperatures for asteroids, the wavelength range that contributes most to this average is between 8 
to 40 microns. This wavelength range contains absorptions featu res due Si-O stretch and bend in silicate mate rials. 
The spectral emissivities are generally between 0.8 and 1.0 (e.g.. 


Salisbury et al 


1991 


Christensen et al. 


20001). As 


the shapes modeled in this work are convex, there is no need to take into account topographic shadowing effect and 
the heating due to the light reflected and emitted by facets on other facets. The effect of roughness on the thermal 
infrared flux is accounted for by adding a spherical-section crater to each surface element of the shape. The crater 
with an opening angle and the crater areal density with respect to the flat part of the surface element pc can be 
va ried from 0 to 90 ° and fr om 0 to 1, respectiv ely, to cover different values of the roughness. Following the procedure 


of 


Delbo’ & Tangal (l2009h and 


Muellen (12012h . we calculate the TPM only for a set of ten roughness models, whose 


parameters are given in Tab. [3] The correspondence between the Hapke’s mean surface slope 0 and the adopted values 
of yc and pc is also given in Tab. [3] Note the degeneracy in the 0 parameter - different combinations of y^ and pc give 
the same 0. 
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Table 3: Ten different values of surface roughness used in the TPM. The table gives the opening angle yc, the crater areal density pc, the Hapke’s 
mean surface slope 6, and our designation. 


Tc 

Pc 

e 

Designation 

0 

0.0 

0.0 

No roughness 

30 

0.3 

3.9 

Low roughness 

40 

0.7 

12.6 

Medium roughness 

40 

0.9 

16.1 

Medium roughness 

50 

0.5 

12.0 

Medium roughness 

60 

0.9 

26.7 

High roughness 

70 

0.7 

27.3 

High roughness 

90 

0.5 

38.8 

High roughness 

90 

0.7 

48.4 

Extreme roughness 

90 

0.9 

55.4 

Extreme roughness 


Shadowing of crater facets on othe r crater facets and mutual heating is taken into account. However, contrary to 
the approach of lPelbo’ & Tangal(l2009h . heat d iffusion is not ex plicitly calculated within craters except for (1865) Cer¬ 


berus. Instead, the analytical approximation of lLagerrosl(ll998h . valid for small solar phase angles, is used. The phase 
angles for some NBAs can reach i n the case of W ISE observations even values of ~90°. The Lagerros approximation 


cannot be used on the night side (ILagerros 


19981) . which, however, corresponds to the half of the surface if a ~9Q°. 


Asteroid (1865) Cerberus is the only one from here studied objects affected by the high phase angle observations, 
thus we have to explicitly calculate the heat diffusion in the craters in this case. For the remaining asteroids, fluxes 
computed by both approaches differ by less than 1%, which is less than are the uncertainties of the observed fluxes 
we use. We favor the Lagerros approximation because of a considerable reduction of the computational time (factor 
of ~40). 

Given the asteroid convex shape and its rotational state, we run the TPM model for different values of the thermal 
inertia T 6 (0,2500) J m“^ s K * and ten combinations of surface roughness (Tab.O. For each value of the surface 
roughness, we run the TPM for the thermal inertia T = 2500 J m“^ s and the Bond albedo A - 0.08, and get 

the first size estimate D. Prior each following TPM run (while k eeping the same surface r oughness), we first compute 


the new value of the Bond albedo A from the equation (see, e.g.. 


Harris & Lagerrosll2002h 


Z)(km) = — lO-O'^", 
Wv 


( 2 ) 


where we use diameter D determined in the previous TPM run, and where t he visible geometric albedo py can 


be expressed via A ^ q py, where q = 0.290 + 0.684G is the phase integral (IBowell et al 


1989h . We adopt the 
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values of absolute magnitudes H and slopes 

Muinonen et al. 

2010; 

Oszkiewicz et al. 

2011) 


201 Ih . Following the procedure of 


Mueller! (l2012h . we then run the TPM 


model with the next (lower) value of F from our grid all the way until F = 0 (always recomputing the A values 
before each following step). We start the iteration with the highest value of F from our grid, because we need 
few steps to reach the realistic combination of A and D. While thermal inertia values for majority of asteroids are 
expected to be significantly lower, sometimes even bellow 10 J m“^ s"^^^ K“', starting the iteration with highest value 
of thermal inertia, could less likely bias the results than starting with the lowest one (i.e., 0 J m“^ s“'^^ K^')- In any 
case, according to our tests for several objects, starting the iteration from zero thermal inertia gives consistent TPM 


results for F > 5 - 10 J m 


-1/2 K- 


Shape models are derived by the lightcurve inversion technique, where a large parameter space is searched; side¬ 
real rotational period, pole orientation, shape, scattering parameters. Due to the uncertainty in the period, the initial 
rotational phase is typically known with an accuracy of +10 degrees. As a result, the orientation of the shape model 
(i.e., the rotational phase <^o) is known to this level of certainty inside the interval of optical lightcurve observations 
used for the lightcurve inversion. When the thermal infrared data are acquired outside this time span, however, the 
uncertainty of 0o propagates proportionally with respect to the time elapsed between the acquisition of the thermal and 
the optical observations. For more details see Appendix B. For most studied asteroids, their expected uncertainties 
of 00 are 6(po < 20°, which implies that it is necessary to include the rotational phase 0o into the TPM optimization 
because even a change of 3^ degrees can significantly improve the fit. We decide to scan 0o with a step of 2° within 
the expected range. 

An example of such sequence of TPM runs for five different values of surface roughness is shown in Fig.lT] where 
we plot the dependence of the reduced values of the TPM fits on the thermal inertia for asteroid (1627) Ivar. 

After scanning the parameter space of the thermal inertia F, surface roughness 0, Bond albedo A and rotational 
phase 00 , we find the solution with the lowest value. The metric is used to find the best-fitting solution takes 
into account directly only the uncertainties of the thermal infrared fluxes. It neglects uncertainties of the shape model 
and the pole orientation. We investigate their importance in thermophysical modeling below. 

We scale the size of the shape model directly by the TPM. After that, we determine the volume equivalent diameter 
D (i.e. the diameter of a sphere with the same volume as the scaled shape model), which is a quantity we present 
throughout this work. 


All-Lagoa et al. 

2014 

Emerv et al. 

2014) 


use the;^^^ statistic to estimate the goodness-of-flt (see, e.g., 
chi-square is computed from the standard chi-square metric defined by Eq. [T]as where the number of 

degrees of freedom v corresponds to the number of data points minus the numbe r of free paramete rs. All solutions 
within < (1 + cr), where cr = V2v/v, are considered indistinguishable (e.g.. jPress et al.lll986h and are used to 


'http://wiki.helsinki.£i/display/PSR/Asteroid+absolute+magnitude+and+slope 


12 






















estimate the uncertainties of the fitted parameters. However, this approach is, in principle, reliable only if;^^^ 


red 


in the (ideal) case of normally distributed and independent errors of the thermal infrared fluxes. Values of ^ 1 
indicate that the model fit has not fully captured the data, but such solutions can still produce realis tic results if care 


Delbo’ & Tanga 


is exe rcised in the analysis. Values of;^^^^^ » 1 indicate a poor fit that should be rejected (e.g., in 
20091 the authors rejected solutions with ~ 8). An alternative is to estimate the fitted parameter uncertainties 


by means of an empirical approach: the error bars include all parameter values corresponding to all solutions with 
^red ^ * (1 + o")- Such Uncertainties often span well the minima in the thermal inertia parameter space for our 

typical values of 1-3, and we always check the appearance of the minima individually to ensure that this is the 


case. 


3.2. TPM with fixed shape models 




100 200 300 400 

Thermal inertia 


500 



Figure 2: Thermal inertia F fits for the best-fitting surface roughness and rotational phase for nine studied asteroids (DAMIT shape models). For 
asteroids (832) Karin, (1472) Muonio and (2606) Odessa, we plot the dependence for both ambiguous pole solutions. The dashed horizontal line 

indicates ;rred = '■ 
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Table 4: Thermophysical characteristics of asteroids derived by the classical TPM with the DAMIT (D in the “pole" 
column) and revised (R) shape models, as well as by the VS-TPM (VS). We provide the asteroid number and name, the 
pole solution, the number of IRAS thermal IR measurements Viras, the number of WISE thermal IR measurements in 
filters W3 Aw 3 and W4 Aw 4 , volume equivalent diameter D, thermal inertia T, visual geometric albedo p\, Hapke’s 
mean su rface slope 6>, rotational phase (^n, reduced chi- square of the best fit absolute magnitude H and slope G 
(AAMS, Muinonen et akll^OlO : Oszkiewicz et al. 2011 ). taxonomical class and average helio centric distance r of the 


asteroids when observed by WISE. Eor the taxonomy, w e primarilly show SMASS II class ( Bus & Binzel 


not available, then the Tholen class (lTholerJl984 


19891) or taxonomy from 


DeMeo & CarrvI (I: 


20021), if 


20lJ. 


Asteroid 

Pole 

Mras 

Aw3 

Nw 4 

D 

[km] 

E 

[Jm-^s-O^K-i] 

pv 

0 

00 

[deg] 

xid 

H 

[mag] 

G 

TAX 

r 

[AU] 

624 

Hektor 

ID 


11 

11 

181!!o 

61 

^■^'^-0.01 

26.7 

354 

3.3 

7.3 

0.33 

D 

5.3 



IR 




1751, 

61 

0.056™ 

26.7 

2 

4.5 







IVS 




186^4 

61 

0.058™ 

26.7 


2.4 





771 

Libera 

ID 

6 

31 

31 

30.5![;^ 

70!g 

0 lA+O.Ol 

27.3 

2 

4.2 

10.3 

0.32 

X 

2.8 



IVS 




321 

oj_3j 

^ *^-0.02 

26.7 


3.4 





832 

Karin 

ID 


12 

12 

16.4!l 

170+90 
^ '^-80 

0 21+0.04 
^■^^-0.01 

38.8 

12 

4.6 

11.1 

0.16 

S 

2.9 



2D 




16.81J 

190+110 

^■^^-0.01 

38.8 

6 

3.5 







IVS 




16.6!l 

qa+150 

^^-90 

A 99+0.04 
^■^^-0.04 

26.7 


1.9 







2VS 




171 

65!1^ 

0 23+0.05 
^■^*^-0.06 

16.1 


1.6 





1036 

Ganymed 

ID 

2 

22 

22 

35.51;^ 

4010 

0 26'*'0 01 
^■^^-0.03 

26.7 

140 

2.8 

9.2 

0.31 

s 

3.9 



IR 




36.01:“ 

351“ 

A 9C9+O.OO8 

26.7 

12 

1.2 







IVS 




-4 

j>2_29 

0 25+0-05 
^ “^-0.03 

26.7 


1.1 





1472 

Muonio 

ID 


10 

9 

8.61:^ 

1+29 

^-1 

0 24+0 05 

3.9 

358 

4.0 

12.4 

0.35 

- 

2.7 


























Table 4: continued. 


Asteroid 

Pole Niras 

Vw3 

Vw4 

D 

r 

Pv 

0 

00 

+L 

H 

G 

TAX 

r 





[km] 

[Jm-2s-02 K-i] 



[deg] 


[mag] 



[AU] 


IVS 



9 1*'' 

25!“ 

A 90+0.09 
^■"^■^-0.06 

16.1 


1.3 






2VS 



9.2*0'3 

0!^° 

A 90+0.09 
^■^■^-0.06 

0.0 


2.2 





1627 Ivar 

ID 

13 

13 

7 44-0.2 
-0.2 

180!™ 

0 26 
^■^^-0.02 

16.1 

4 

1.0 

12.6 

0.33 

S 

2.1 


IR 



7 4+0.1 

ioo!2o 

-0.005 

12.0 

10 

0.8 






IVS 



^■Kl 

100+ 120 

rou 40 

u.z,jj_00i4 

12.0 


0.8 





1865 Cerberus* 

ID 

10 

9 

-1.3 

-1250 

-0.17 

16.1 

350 

23.1 

16.6 

0.37 

s 

1.1 

* 

IR 



-1.2 

-400 

-0.215 

55.4 

352 

16.2 





* 

IVS 



-1.2 

>300 

-0.30 

55.4 


6.3 





1980 Tezcatlipoca 

ID 

17 

15 

5 0** * 

60+390 
uo 5g 

A 23+0.07 
^■^*^-0.08 

16.1 

356 

4.6 

13.6 

0.18 

SI 

2.3 


IR 




^0+440 

^^-42 

0 29+o *)5 

16.1 

350 

3.5 






IVS 



5.4!0.7 

220!280 

^■^^-0.07 

48.4 


2.9 





2606 Odessa 

ID 

9 

9 

17!^ 

100!™ 

A 10+0.06 
^■^^-0.02 

12.0 

348 

8.9 

11.7 

0.21 

X 

3.5 


2D 



16!^ 

qc+85 

^■^^-0.03 

12.6 

348 

5.6 






IVS 



20!^ 

ioo!^f 

^■^^-0.02 

12.0 


1.6 






2VS 



17*2 

90!«° 

0 13+0.04 
^■^■^-0.04 

12.0 


1.5 






*Poor TPM fit 






We take shape models of four near-Earth asteroids (1036) Ganymed, (1627) Ivar, (1865) Cerberus and (1980) Tez- 
catlipoca, of four main-belt asteroids (771) Libera, (832) Karin, (1472) Muonio and (2606) Odessa, and of the Jovian 
Trojan asteroid (624) Hektor from the publicly available DAMIT database. In Tab. [T] we list rotational parameters, 
information about the photometry used for the shape model determination, and the references to the original publica¬ 
tions. These solutions have well defined sidereal rotational period values (the accuracy corresponds to the last decimal 
place in period values in Tab.[T]), and usually two symmetric pole orientations with similar ecliptic latitudes p and dif¬ 
ference in ecliptic longitude d of ~ 180° (the so-called pole ambiguity, a typical example is asteroid (832) Karin). 
However, due to the more various observing geometries of NBAs than of the main-belt asteroids or Jovian Trojans, 
the pole ambiguity is not present in our sample of NBAs. Moreover,!] 


Marchis et al 


(12014h removed the pole ambiguity 


for asteroid (624) Hektor thanks to the disk-resolved images, thus we use only the preferred solution (see Tab[T]i for 
the TPM. 


We run the TPM for each studied asteroid and for various initial values of thermophysical parameters as described 
in Sect. 13.11 In Tab. H] we list thermophysical parameters of the best-fitting solutions -T, D, 6 and the geometric 
visible albedo pv (instead of the Bond albedo A). 

We obtain best-fitting TPM solution with;^f^^^ ~ 1 only for asteroid (1627) Ivar. This solution can be considered 
formally acceptable. 

Bor asteroids (624) Hektor, (771) Libera, (832) Karin, (1036) Ganymed, (1472) Muonio, (1980) Tezcatlipoca and 
(2606) Odessa, we obtain TPM fits with higher values in the range of 3-9, which means that our model is not fully 
reproducing the thermal observations. The extreme case is asteroid (1865) Cerberus with^^^^ > 20, which means that 
we failed to constrain any of its properties. 

The best-fitting thermophysical properties and their uncertainty estimates are included in Tab. |4| 

The convergence of the solution in the thermal inertia parameter space for all studied asteroids is shown in 
Pigs. |2^,b,c,d, where the curves correspond to the best-fitting surface roughness and rotational phase for each 
value thermal inertia. We show both pole solutions when they are available. This figure shows that there is a more 
or less prominent minimum for most asteroids, and therefore the thermal inertia is usually well constrained 
(e.g., (1036) Ganymed, (1472) Muonio or (2606) Odessa). The thermal inertia for asteroid (624) Hektor is the best 
constrained one (T < 30 Jm^^s '^^K '), but in general, the range of acceptable thermal inertia values is broader 
(see Tab. nil. Bor asteroid (1980) Tezcatlipoca, the thermal inertia is poorly constrained providing only an upper limit 
of <500Jm-2s-O2K-i. 

The highx^^^ values could be, in principle, explained by: (i) shape model uncertainties (i.e., convex shape, pole 
orientation) dominate over the flux uncertainties, (ii) systematic uncertainties due to various assumptions of our TPM 
model are biasing modeled fluxes (for example, a finite grid of searched parameters, or the usage of a specific rough¬ 
ness model), or Anally (iii) an underestimation of the uncertainties of the WISE thermal infrared measurements (in¬ 
cluding systematics in calibration and outliers). These are discussed in the following subsections. In SectjH we present 
and apply the varied-shape TPM method, which shows that accounting for the uncertainties in the shape model and 
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the pole orientation is essential for the determination of reliable thermophysical properties by TPM. 

TPM modeling of asteroid (1036) Ganymed is unique because the rotational phase (^o is not constrained at all. 
The shape model is based on photometric data from only one apparition in 1985, so the uncertainty of the sidereal 
rotational period of 0.001 hour for the WISE observations (year 2010) corresponds to a rotational uncertainty up to 
~ 540 degrees. 

In all three cases with ambiguous pole solutions, the TPM produce different fits (more than one in the best-fitting 
Tred values), which suggests that one pole solution is preferred, however, the thermal inertia values are consistent. 

The fits to the thermal infrared data for selected/all asteroids are shown in the Supplementary material. 


3.3. Revised shape models 



Thermal inertia [Jm 



Figure 3: Thermal inertia F fits for the best-fitting surface roughness and rotational phase for DAMIT (blue) and revised (red) shape models of 
asteroids (624) Hektor, (1036) Ganymed, (1627) Ivai; (1865) Cerbems and (1980) Tezcatlipoca. The dashed horizontal line indicates = 1. A 
color version of the figure is available in the electronic version of the journal. 


For some asteroids new optical lightcurves were acquired after their shape models were determined and included 
in the DAMIT. Thus we have the opportunity to revise the shape models of asteroids (624) Hektor, (1036) Ganymed, 
(1 627) Ivar, (1 8 65) Cerberus and (1980) Tezcatlipoca. All revised shape models, derived following the procedure 


of 


Hanus et al 


(12011 


2013bl) . are based on combined dense- and sparse-in-time photometric data sets. Rotational 


parameters and information about the photometry are included in Tab. [T] 

The typical uncertainty of the pole orientation, which depends on the amount and quality of available photometric 
data, is in the ecliptic coordinate frame (5-10°)/cosy6 in longitude A and 5-10° in latitude [3 (see also Fig.|4] where we 
map a typical pole uncertainty for all studied asteroids). The same applies to the DAMIT models we already used in 
Sect.im 

We run the TPM with the revised shape models as fixed inputs as in the previous section and we present the derived 
thermophysical parameters in Tab. |4] The best fitting values of TPM fits do not change significantly but the quality 
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of the fits is generally improved (see the comparison in Figs. [3^,b). The only exception is (624) Hektor, for which 
the TPM fit with the revised shape model is slightly worse, but consistent in the best-fitting values of thermophysical 
properties. 

These results indicate that the introduction of revised shape models affects the goodness-of-fit. For asteroid 
(1627) Ivar, the TPM with the revised shape model reduced the uncertainty of the thermal inertia value by a fac¬ 
tor of ~3. In turn, we do not significantly lower the values for asteroids (624) Hektor, (1865) Cerberus and 
(1980) Tezcatlipoca, thus we still do not ideally reproduce the thermal infrared data. On the other hand, the TPM fit 
of asteroid (1036) Ganymed has now best-fitting value close to one. 

One question that still remains is whether some improvements in the fitting technique could be the main cause of 
omx^ed values being higher than one. In the next subsection we investigate this and other related issues. 

3.4. Additional model uncertainties 

Because we run TPM with thermophysical parameters (thermal inertia F, surface roughness 0, rotational phase 0o) 
from a finite grid of values, our best-fitting solution could not correspond to a real global minimum in the searched 
parameter space; however, it should be close. By using a finer grid of parameter values, we should be in principle 
able to find a solution with a lower By varying F and <po, we are able to reduce the;^f^ for all studied asteroids 
typically by only few percent, which is a rather marginal improvement. This means that our grid in these parameters 
is sufficiently fine. 

To deal with the surface roughness, we use the Lagerros analytical crater approximation, which is computationally 
less demanding than performing a full heat diffusion computation within the craters. The differences in the fluxes 
produced by these two models are usually considerably smaller than 1 %, which are always lower than the uncertainties 
of the observed fluxes by a factor of at least five. As a result, the quality of the fit should be affected only marginally. 

On the other hand, we use only ten different values for the surface roughness 9 in the modeling (see Tab.O. In the 
practical application, we rather use the opening angle yc and the crater areal density pc than the Hapke’s mean slope 
0. By using more combination of and pc, we are able to reduce thsx^ values of the TPM fits by few percent, which 
is again a marginal improvement of the fit. 

We also do not account for the uncertainty in the input H and G values. According to our tests (TPM with different 
H values), a change of +0.5 mag in H is compensated by the change in Bond albedo A. However, the size remains 
similar (see Eq. |2l), as well as the best-fitting values of thermal inertia, surface roughness, or indeed the resulting 
minimum 

Our results call for additional investigation of the role of the shape model uncertainties in the TPM, which we 
perform by using our novel VS-TPM method in the next section. 
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Figure 4: Orientations of the spin axes of all varied shape models for here studied asteroids. A color version of the figure is available in the 
electronic version of the journal. 


19 




Figure 5: Six shape models of asteroid (1627) Ivar derived by the lightcurve inversion method by bootstrapping the optical photometric data sets. 
First two images are equator-on views rotated by 90° and the third one is a pole-on view. 


4. Varied shape TPM 

4.1. VS-TPM - method 


It is known that the disk-integrated photometry is not significantly influenced by the non-convexities and shape 
features, and that the conv ex inversion is usually s table and produce similar results even if we reasonably vary the 


photometric data (see e.g., 


Kaasalainen et al 


200 ih . As a result, most of the available shape models are convex. In 


the case of very precise thermal measurements (such as from the WISE satellite), the shape features together with 
the pole uncertainty could play a more significant role for the reliable determination of thermophysical properties, 
because they could considerably influence the temperature profile on the surface and thus the thermal infrared fluxes. 
To investigate this behavior, we introduce here a procedure that varies both the shape model and the pole orientation, 
while keeping the shape convex. Each varied shape model can be then used as input for the TPM modeling. We call 
this method varied shape TPM (VS-TPM). 

Instead of varying the shape models directly^, we bootstrap the photometric data sets and use the convex inversion 
method to derive varied shape models. The advantage of this approach is that it is based only on the data and it also 
maps uncertainties in the pole direction. 

The whole procedure of investigating the stability of the TPM on the shape model variations can be divided into 
the following steps: 


*For example by deforming the polyhedron, by introducing concavities, or by more sophisticated methods such as transforming the polyhedron 
into the Gaussian image, varying the facets and transforming it back, or by representing the shape by spherical harmonic series and varying the 
coefficients. 
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1. Bootstrapping of the photometric data. We randomly select a number of dense lightcurves from the original pho¬ 
tometric data set equal to the number it contains, i.e., if the original data set contains ten individual lightcurves, 
we randomly select ten lightcurves from that data set. As a result, we can have some lightcurve multiple times 
in the new data set, and some of the lightcurves can be missing. This, in principle, corresponds to a different 
weighing of the lightcurves. For the sparse data, the modification procedure is similar: we randomly choose 
individual measurements from the pool of original data, while keeping the same number of observations in 
the whole sparse-in-time lightcurve. D ata from different astrometric sources are treated separately. Generally 


speaking, we use the bootstrap method (IPress et al.lll986h . 


2. Shape model determination. We use the randomly selected photometric data set to derive the varied shape model 
by the lightcurve inversion. The original rotational state solution serves as an initial guess for the optimized 
parameters. The rotational state is usually close to the original one as can be seen in Fig. |4] where we show 
the orientations of the spin axes of varied shape models of all nine asteroids studied here. The difference in the 
rotational state corresponds to the expected uncertainty, which is typically 5-10° in the pole orientation. The 
shape appearance is similar, as can be seen in Fig. [S] where we show several varied shape models of asteroid 
(1627) Ivar based on modified photometric data sets. 


3. TPM modeling. We perform the TPM optimization scheme for each varied shape and its rotational state the 
same way as for the original shape model. 


4. We repeat steps 1-3 to obtain a desired statistical sample. We run VS-TPM for each pole solution individually. 


All the varied shape models fit the visible lightcurves equally well. Usually, the dense lightcurves are not weighed 
when the shape models are computed by the lightcurve inversion, but it is obvious that individual lightcurves do not 
have equal quality. On the other hand, the sparse data set is weighed with respect to the dense data to penalize its 
lower photometric accuracy. Our method not only naturally varies the lightcurve weights, but also selects subsets of 
the photometric data by ignoring some lightcurves, and thus produces sets of slightly different shape models whose 
rotational parameters sample the underlying uncertainties of the optical data. 


4.2. VS-TPM - application to nine asteroids 

In Sects. and 13.31 we already illustrate for several asteroids that the accuracy of the shape and its rotational 
state should be considered in the TPM modeling. In particular, in the case of (1865) Cerberus, the classical TPM 
approach leads to slightly different fits when the revised or the DAMIT shapes are used, indicating that the solution 
does not appear stable against variations of the shape and the pole orientation. Additionally the high values of 
the TPM fits for asteroids (1472) Muonio, (1865) Cerberus, (1980) Tezcatlipoca and (2606) Odessa indicate that the 
current shape models are not able to reproduce carefully the observed thermal fluxes. We thus apply our novel method 
VS-TPM to investigate the stability of the TPM solution against the shape and pole uncertainties and to estimate the 
contribution of the shape model uncertainty to the;^^^^^ for all nine studied asteroids. 
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Figure 6: Distribution of the best-fitted thermai inertia vatues (red futt circtes) determined for the varied shape modets of the Jovian Trojan 
(624) Hektor and main-beit asteroids (771) Libera, (832) Karin, (1472) Muonio and (2606) Odessa, together with the convergence in the thermal 
inertia (black lines). The second ambiguous pole solution is represented by blue asterisks. We also show the thermal inertia solutions with the 
corresponding error bars derived by the classical TPM for both revised (blue) and the DAMIT (black) shape models. The dashed horizontal line 
indicates = 1- Note the ditferent scales for the thermal inertia. A color version of the figure is available in the electronic version of the journal. 
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Figure 7: Distribution of the best-fitted thermal inertia values (red full circles) determined for the varied shape models of four near-Earth asteroids 
(1036) Ganymed, (1627) Ivar, (1865) Cerbenis and (1980) Tezcatlipoca, together with the convergence in the thermal inertia (black lines). We 
also show the thermal inertia solutions with the corresponding error bars derived by the classical TPM for both the revised (blue) and the DAMIT 
(black) shape models. The dashed horizontal line indicates = i- Note the different scales for the thermal inertia. A color version of the figure 
is available in the electronic version of the journal. 
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We create for each studied asteroid a sample of 29 varied shape models and their corresponding rotational states 
following the scheme described in Sect. 14. II Together with the original shape model, we have a statistical sample of 30 
slightly different shape models for each asteroid. In Fig. 01 we show the distribution of spin vectors of the varied shape 
models for all nine studied asteroids. They are consistent within their expected uncertainties. The typical dispersion 
is -10-20°. 


Thereafter, we run the TPM for each varied shape model and pole solution and show the resulting thermophysical 
parameters in Tab01 In Figs. |6] and I?) we plot the dependence of thex^e^ versus thermal inertias (gray dashed lines) 
for all varied shape models and highlight each minimum; red full circles or blue asterisks, corresponding to different 
ambiguous pole solutions. We also include the solutions from Sect. 13.21 with their uncertainties for both revised (if 
available) and DAMIT shape models. 

The x^g^ of the solutions of the varied shape TPM reach values ~ 1 in five out of nine cases, which shows that we 
are reliably fitting the thermal infrared data and that the VS-TPM supersedes the classical approach. This is because 
varied shape models, even though they are indistinguishable in terms of reproducing the optical photometry, are clearly 
different in the thermal infrared. Note that this method is not an optimization of the shape, rather it is a way to map the 
uncertainties in the TPM analyses related to the uncertainties of the shape and the pole orientation given by the optical 
photometry. A technique that truly optimizes the shape (still as convex), the rotational state and the thermophysical 
properties simultaneo usly from optical an d the thermal infrared photometry is under development (multiple source 


data inversion method. 


Durech et alJl2014t). 


We observe various behaviors for the distribution of the best TPM solutions in our sample of asteroids. In the cases 
of asteroids (624) Hektor, (1036) Ganymed and (1627) Ivar, the chi-square curves and their corresponding minima 
present qualitatively similar behavior and lead to statistically indistinguishable solutions. This suggests that the TPM 
results are not strongly dependent on the varied shape models, likely because the optical photometric data sets for 
these asteroids are particularly large and, therefore, their shapes are well constrained. Not surprisingly, the range of 
acceptable thermal inertia values is larger than suggested by the classical TPM analysis. 

The most intriguing result is the fact that we also obtain TPM fits with a large range of x^^,^ values for most 
of the studied asteroids. This suggests that TPM modeling with individual varied shape models differ considerably, 
and sometimes the best fitting solutions have x^^^ values that differ by a factor of 5 while keeping thermal inertia 
determinations consistent, as can be seen in Fig.|6]for asteroids (771) Libera or (1472) Muonio. The most extreme 
case is asteroid (2606) Odessa, with;^'^^^ between 2 and 12. Contrary to the classical TPM approach, the application of 
the VS-TPM has enabled us to obtain satisfactory fits (with;\f^^^ ~ 1) for asteroids (1472) Muonio and (2606) Odessa. 

On the other hand, the best thermal inertia solutions within the varied shape models of asteroid (832) Karin range 
from zero to two hundred, some of which are not consistent with the values determined by the classical approach. 
This is a warning that fixing a single shape model with the classical TPM and neglecting the uncertainties in the shape 
and the pole orientation can strongly bias the solution. In turn, the VS-TPM cannot fit the thermal infrared data of 
asteroid (1865) Cerberus. The distribution of the best fitting thermal inertia solutions (red full circles in Fig.lT]) as well 
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as their correspendingfor (1865) Cerberus span the largest range of values. This shows a very strong dependence 
on the shape model and rotational state, probably enhanced by the high phase angle of the observations (a ~ 90) and 
the extreme elongation of the shape. 

To sum up, derived thermal inertia values of asteroids (624) Hektor, (771) Libera, (1036) Ganymed, (1472) Muo- 
nio, (1627) Ivar and (2606) Odessa seem to be well constrained, contrary to those of (832) Karin and (1980) Tez- 
catlipoca. The TPM solution for (1865) Cerberus is still poor. Contrary to the classical TPM, the VS-TPM approach 
does not always help reject one of the ambiguous pole solutions. In particular, there is not a preferred solution for 
asteroids (832) Karin, (1472) Muonio and (2606) Odessa. 

To estimate thermophysical properties and their corresponding uncertainties based on the VS-TPM, we consider 
the confidence interval given by the 21 (~ 68 %) TPM solutions with the lowest best-fittingThis criterion allows 
us to exclude the most extreme bootstrapped shapes from the analysis, and, simultaneously, to account for the TPM 
parameter ranges spanned by the varied shape models. We calculate the uncertainties of the fitted parameters by 
the empirical approach the same way as in Sect. 13.21 for each single varied shape: xfed ^ * (1 + V^/i^)- By 

finding the maximum and minimum values allowed by the individual uncertainties within these 21 cases, we obtain 
the final uncertainties quoted in Tab.|4] Note that the TPM fit for (1865) Cerberus is poor, thus we do not report the 
uncertainties; however, we are still able to estimate a lower limit of the thermal inertia. 


4.3. Individual cases 


Mueller et al. 


(624) Hektor. The low thermal inertia of this D-type asteroids is no t surprising, it is sli ghtly smaller than the one 
of another Trojan asteroid (617) Patroclus (20+15 Jm“^ s“'^^ K“', 

(62+37 Jm-2 s-i/^K 
for (3063) Makhaon, 


for (1173) 

Anchises, 7+7 Jm ^ 

S-l/2]< 

Horner et al. 

2012 

Fernandez et al. 

2003 ) 


20101) . or other Jovian Trojans 

ac Ti rr,-2 c-1/2 ir-l 


Note that Hektor is a binary minor body consisting of a primary component and a moonlet with a size ratio of 
XXX. Thus, infrared fluxes are dominated by the thermal emission of the primary and the moonlet can be neglected 
in the TPM analysis. 


Our volume equivalent diameter of D = 186^^34 km is Icr larger than the diameter reported by 


Gray et al. 


(l 2012 bl) 


based on the NEATM, but lower th an the diameter base d on the AKARI data (230 km, based on a standard thermal 
model. 


Usui et al 


201 Ih . Moreover. 


Marchis et al 


(12014ft used their own convex shape model to fit the contour profiles 


of adaptive optics images acquired by Keck telescope and derived a volume-equivalent diameter of D — 270 + 22 km. 
We also perform the TPM analysis of the W ISE data with t he con vex shape model of y 
and still obtain a significantly smaller size. 


Marchis et al 


Marchis et al 


(l2014ft as input 


(12014ft give some evidence that the primary component 


of (624) Hektor has a bi-lobed shape. If this is the case, as the convex shape model usually corresponds to the 
convex hull of the real shape, the diameter estimated from TPM could be, in principle, both over- or underestimated 
depending on the observing/viewing geometry. In our particular case, the observing configuration during the (almost) 
maximum elongation is such that the projected area of the convex model would be bigger and flatter than that of the 
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bi-lobed shape model. Given the low phase angle of the observations the model temperatures of the convex shape 
are expected to be higher than those of a non convex shape. Consequently, the best fitting diameter becomes smaller 
to compensate for the higher temperatures on the visible hemisphere. While we caution that more accurate adaptive 
optics images, and/or occultation data are likely necessary to investigate the shape of this body, we point out that the 


value of pv = 0.058 /q™ corres ponding to D 


trojans of 0.07+0.03 reported by 


Gray et al 


190 km is in very good agreement with the average albedo of large 


(1201 ih . On the other hand, the value of pv = 0.029, corresponding to 


D ~ 270 km is significantly lower. 

Jovian Trojan (624) Hektor is another asteroid with an emiss ion band at ~ 10 pm and low thermal inertia (such as 
(24) Themis, (1) Ceres, or (617) Patroclus, Vernazza et aklLoi^ ). Such low thermal inertia could be consistent with a 


very fine and porous surface material (fluffy material). 

(1036) Ganymed. The photometric data set for this asteroid is exceptionally large and the TPM produ ces thermal 
inertia of T = 35 / 29 , which is comparable to the value for the Moon’s fine regolith (ISnencer et alJll989r) . Such low 


value of thermal inertia is not common among NBAs (indeed, Ganymed is a NBAs with t he second lowest B so far), 

. iDelbo’ et alJl2007 ). 


for which most determinations are of the order of some hundreds J m ^ s K * (see, e.g.. 


However, we note that the thermal conductivity in the regolith is temperature dependent (IKeihm 


19841) . and so 


is thermal inertia. Under the assumption that heat is transported in the regolith mainly by radiative conduction be- 


Kuhrt & Giese 

19$9 

Jakoskv 

1986) 


Jakoskvlll986l) . In this case B oc and, because T oc then B oc r where r is the helio¬ 


centric distance of the body. As Ganymed was observed at r ~ 3.7 AU we expect that its thermal inertia value at 1 AU 
would be around three times higher, namely ~100-150 Jm“^ s"'^^ K^'. Accurate modeling of the heat transfer in the 

- ir^ 

regolith (IGundlach & Blumll2013l) shows that the heat conduction between touching regolith grains is also important 
and the temperature depe ndence of the conductivity might depart from the pure radiative term. Using the model of 
Gundlach & BlumI (120 1 ■3h and assuming S-type like thermal properties we can estimate that the thermal inertia can 
increase by a factor between 1.3 and 2.7 as Ganymed moves from 3.7 to 1 AU depending on the average grain size 
and regolith packing fraction. We conclude that the thermal inertia value at 1 AU of Ganymed is between 45 and 
100 J m“^ s“'^^ which is still am ong the lowest valu es measured for NBAs (the lowest value so far is B = 36^20 
J vaT^ s K ' for asteroid 1950 DA, 


Rozitis et al.ll2014l). 


Another explanation for the low thermal inertia could be the unusually large size of this NBA: D - km. This 
is significantly larger than the other NBAs with thermal inertia determin ations (usually of th e order of 1 km). This is 


also consistent with the size vs. thermal inertia dependence proposed by 


Delbo’ et al 


(l2007l) . 


The lower-than-average value for the thermal inertia of Ganymed is likely due to the presence of a finer than 
average surface regolith, implying that this asteroid has a surface that has been exposed longer than the average 
NBA to the micrometeorite bomba rdment, which is one of the mechanisms that has been claimed to be capable 


of comminuting rocks on asteroids (IHorz & Cintala 


19971) . A more convincing explanation is thermal cracking: it 
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has been recently shown that thermal fragmentation induced by the diurnal temperatur e variations br eaks up rocks 


larger than a few centimeters more quickly than do micrometeoroid impacts on NBAs (IDelbo’ 


20141) . Because the 


process of thermal fragmentation is strongly dependent on the value of of the diurnal temperature excursion, which 
increases with decreasing perihelion distance, we expect NBAs with low perihelion distances to have finer regolith and 


NBA s according to 


Marchi et al. 

o 

o 

) (see their Big 4). Analysis of the orbital history of Ganymed ( 

Marchi et al. 


2009h shows that this asteroid has 98% probability to have had a perihelion distance smaller or equal to 0.1 AU (!) in 


the past for at least 4,000 years, which implies that thermal cracking was very efficient for this asteroid. 

The geometric visible albedo is pv = 0.25;|^gg3, which is a typical value for an S-type asteroid. Based on the value 
of B, we expect that the regolith grain sizes are comparable to those on the Moon. 


5. Discussion 

Asteroids shape models are determined from optical lightcurves that usually cover only several apparitions, and 
therefore a limited number of observing geometries. The shape model then predicts well the lightcurves for the 
apparitions covered by the observations but could be less accurate for other apparitions, especially when the geometry 
of observations is significantly different. In such cases, parts of the surface that are exposed to the observer could not 
be realistic. If the shape model is not based on dense-in-time photometry from the same apparition as the WISB data 
were observed in, we should expect that the computed fluxes are less accurate/realistic. 

Most TPM works so far have relied on convex shape models for several reasons: (i) the majority of shape model 
determinations are convex, (ii) the high quality optical photometric data (i.e., with high photometric accuracy and/or 
acquired at high phase angles) necessary for a non-convex shape model determination is available for only very few 
asteroids, moreover, there are usually no WISB thermal infrared data for them, because filters W3 and W4 saturate 
for big asteroids, (iii) radar observations, which are typically used for non-convex shape model determinations, are 
limited to only few largest and closest MBAs or close flybys of small NBAs, or (iv) the convex shape model usually 
fits the optical data to the level of the noise, thus optimization with a non-convex model does not provide meaningful 
results because the disk-integrated photometry usually contains little information about the non-convexities. 

On the other hand, the convex approximation could be suffici ent because convex and non-convex models of aster¬ 


oid (1620) Geographos give consistent thermophysical sol utions (iRozitis & Gre^l2014l) . and thermophysical proper¬ 


ties of asteroid (21) Buteti a based on a convex s hape model (ICarrv et al 
the Rosetta/VIRTIS data ( Coradini et alJl^OH ). It was also shown in fMarchis et al 


2012 ) are consistent with the findings base d on 


(l2006l) : 


Hanus et al 


(l2013al) that 


the disk-resolved images acquired by the 8-lOm class telescopes equipped with AO systems usually well correspond 
to the asteroid’s 2D shape projections, thus the shape models globally well represent the real asteroid’s appearance. 
The fact, that in some cases we are not able to reproduce well the observed fluxes by the TPM with such convex shape 
models (i.e., with those that are reliable global representations of the true shapes), or that individual varied shapes 
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give significantly different values of the TPM fits, provides the evidence that the observed thermal emission is 


indeed more sensitive than optical observations to topographic features and concavities that do not significantly alter 
the global shape. 

Our results show that the shape model (together with the pole orientation) is a significant limiting factor for the 
goodness of the TPM fit. We expect this effect to be more prominent when the shape model is based mainly on 
sparse-in-time opticaL data. We also show the stability of the values of the physical parameters derived from TPM 
with respect to changes to the shape model and the pole orientation. This raises our confidence that thermophysical 
properties are in general reliable. 

The VS-TPM also shows that the quoted errors based on the classical TPM approach, which account only for 
statistical uncertainties of the observed infrared fluxes, should be interpreted as minimum estimates of the real uncer¬ 
tainties. The VS-TPM, which accounts for the uncertainties in the shape and pole orientation of the asteroid, is a very 
viable method to estimate realistic error bars of the parameter values. However, at the moment and for the foreseeable 
future, our shapes are and will be based on optical lightcurve inversion, and thus remain convex. Probably,values 
would be lowered if it was possible to derive a non-convex shape model from lightcurve inversion. 

Finally, a model capable of optimizing photometric and thermal infrared data simultaneously (multi-data inversion 
technique) should be a promising step forward to wards the improvem ent of (still convex) asteroid shapes using thermal 


infrared data together with optical observations (iDurech et alJl2014l) . 


6. Conclusions 

We present a novel method to investigate the importance of the shape model and the pole orientation uncertainties 
in the thermophysical modeling - the varied shape TPM (VS-TPM). 

We apply the VS-TPM to nine asteroids and reveal the strong dependence of the TPM fit on the shape model 
and the pole orientation uncertainty for several asteroids. The best-fitting parameters are presented in Tab. |4] and 
discussed in Sect. 14.21 From this table, one can see that the uncertainties of these properties derived by the classical 
TPM method are usually underestimated. In most cases, the best-fitting values of the physical parameter derived from 
the TPM analysis are consistent between the classical TPM and the VS-TPM. However, there are exemples where 
VS-TPM shows that neglecting the uncertainties in the shape model and the rotational state can result in biased values 
of these parameters (e.g., the thermal inertia for the asteroid (832) Karin). 

Furthermore, the VS-TPM allows us to find a significantly better TPM solution than the classical approach (e.g, 
(1472) Muonio and (2606) Odessa). 

Additionally, the obvious differences between the TPM fits found for some of the ambiguous pole solutions based 
on the standard approach are not that prominent after applying the VS-TPM. This suggests that one should be very 
cautious when making conclusions based on the classical results. 
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Based on our findings, we recommend to always consider the uncertainties of the shape model and its pole orien¬ 
tation in the thermophysical modeling. 

With the tools developed in this work, we are now ready to exploit the WISE catalog and determine thermophysical 
properties, especially thermal inertias, for hundreds of asteroids with convex shape models. 

Acknowledgements 

The work of JH, MD and VAL was carried under the contract 11-BS56-008 (SHOCKS) of the French Agence 
National de la Recherche (ANR), and JD has been supported by grants GACR P209/10/0537 and 15-04816S of the 
Czech Science Foundation. VAF acknowledges support from the project AYA2011-29489-C03-02 (MFC, former 
Spanish Ministry of Education and Science). The computations have been done on the “Mesocentre” computers, 
hosted by the Observatoire de la Cote d’Azur, and on the computational cluster Tiger at the Astronomical Institute of 
Charles University in Prague (http: //sirrah, troja.mff. cuni . cz/tiger). 

This publication uses data products from NEOWISE, a project of the Jet Propulsion Faboratory/California Institute 
of Technology, funded by the Planetary Science Division of the NASA. We made use of the NASA/IPAC Infrared Sci¬ 
ence Archive, which is operated by the Jet Propulsion Faboratory, California Institute of Technology, under contract 
with the NASA. 

References 

References 

All-Lagoa, V., Lionni, L., Delbo, M., et al. 2014, Astronomy and Astrophysics, 561, A45 
Binzel, R. R 1987, Icaiiis, 72, 135 

Bottke, J. W. R, Vokrouhlicky, D., Rubincam, D. R, & Nesvomy, D. 2006, Annual Review of Earth and Rlanetary Sciences, 34, 157 
Bowell, E., Hapke, B., Domingue, D., et al. 1989, in Asteroids II, ed. R. R. Binzel, T. Gehrels, & M. S. Matthews, 524-556 
Bus, S. J. & Binzel, R. R. 2002, Icarus, 158, 146 

CaiTy, B., Kaasalainen, M., Merline, W. J., et al. 2012, Rlanetary and Space Science, 66, 200 
Chernova, G. R, Kiselev, N. N., Krugley, Y. N., et al. 1995, Astronomical Journal, 110, 1875 
Christensen, R. R., Bandfield, J. L., Hamilton, V. E., et al. 2000, Journal of Geophysical Reseai'ch, 105, 9735 
Coradini, A., Capaccioni, R, Erard, S., et al. 2011, Science, 334, 492 

Dahlgren, M., Lagerkvist, C.-L, Fitzsimmons, A., & Williams, I. R. 1991, Monthly Notices of the Royal Astronomical Society, 250, 115 
Delbo’, M. 2004, RhD thesis - Freie Univesitaet Berlin, 1 
Delbo’, M. 2014, Nature 

Delbo’, M., deirOro, A., Harris, A. W., Mottola, S., & Mueller, M. 2007, Icams, 190, 236 
Delbo’, M. & Tanga, R. 2009, Rlanetary and Space Science, 57, 259 
DeMeo, R E. & Carry, B. 2013, Icarus, 226, 723 

Detal, A., Hainaut, O., Rospieszalska-Surdej, A., et al. 1994, Astronomy and Astrophysics, 281, 269 
Dunlap, J. L. & Gehrels, T. 1969, Astronomical Journal, 74, 796 


29 


Durech, J., Hanus, J., Delbo, M., Ali-Lagoa, V., & Carry, B. 2014, in AAS/Division for Planetaiy Sciences Meeting Abstracts, Vol. 46, 
AAS/Division for Planetary Sciences Meeting Abstracts, #509.11 
Durech, J., Sidorin, V., & Kaasalainen, M. 2010, Astronomy and Astrophysics, 513, A46 
Durech, J., Voki'ouhlicky, D., Baransky, A. R., et al. 2012, Astronomy and Astrophysics, 547, AlO 
Emery, J. R, Fernandez, Y. R., Kelley, M. S. R, et al. 2014, Icarus, 234, 17 
Emery, J. R, Sprague, A. L., Witteborn, F. C., et al. 1998, Icaius, 136, 104 
Fernandez, Y. R., Sheppard, S. S., & Jewitt, D. C. 2003, Astronomical Journal, 126, 1563 
Grav, T., Mainzer, A. K., Bauer, J., et al. 2012a, Astrophysical Journal, 744, 197 
Grav, T., Mainzer, A. K., Bauer, J., et al. 2011, Astrophysical Journal, 742, 40 

Grav, T., Mainzer, A. K., Bauer, J. M., Masiero, J. R., & Nugent, C. R. 2012b, Astrophysical Journal, 759, 49 

Gundlach, B. & Blum, J. 2013, Icarus, 223, 479 

Hahn, G., Magnusson, R, Harris, A. W., et al. 1989, Icarus, 78, 363 

Hainaut-Rouelle, M.-C., Hainaut, O. R., & Detal, A. 1995, Astronomy and Astrophysics Supplement Series, 112, 125 

Hanus, J., Durech, J., Broz, M., et al. 2011, Astronomy and Astrophysics, 530, A134 

Hanus, J., Marchis, E, & Durech, J. 2013a, Icaims, 226, 1045 

Hanus, J., Durech, J., Broz, M., et al. 2013b, Astronomy and Astrophysics, 551, A67 

Hapke, B. 1984, Icarus, 59, 41 

HaiTis, A. W. 1998, Icarus, 131, 291 

Harris, A. W. & Lagerros, J. S. V. 2002, Asteroids III, 205 

HaiTis, A. W. & Young, J. W. 1989, Icarus, 81, 314 

Hartmann, W. K. & Cruikshank, D. R 1978, Icarus, 36, 353 

Higgins, D., Pravec, R, Kusnirak, R, et al. 2008, Minor Planet Bulletin, 35, 123 

Hoffmann, M. & Geyer, E. H. 1990, Acta Astronomica, 40, 389 

Horner, J., Muller, T. G., & Lykawka, P. S. 2012, Monthly Notices of the Royal Astronomical Society, 423, 2587 

Horz, F. & Cintala, M. 1997, Meteoritics and Planetary Science, 32, 179 

Ito, T. & Yoshida, F. 2007, Publications of the Astronomical Society of Japan, 59, 269 

Jakosky, B. M. 1986, Icarus, 66, 117 

Jarrett, T. H., Cohen, M., Masci, R, et al. 2011, Astrophysical Journal, 735, 112 

Kaasalainen, M., Pravec, R, Ki'ugly, Y. N., et al. 2004, Icarus, 167, 178 

Kaasalainen, M. & Torppa, J. 2001, Icarus, 153, 24 

Kaasalainen, M., Torppa, J., & Muinonen, K. 2001, Icarus, 153, 37 

Kaasalainen, M., Torppa, J., & Piironen, J. 2002a, Astronomy and Astrophysics, 383, L19 

Kaasalainen, M., Torppa, J., & Piironen, J. 2002b, Icarus, 159, 369 

Keihm, S. J. 1984, Icarus, 60, 568 

Kuhrt, E. & Giese, B. 1989, Icarus, 81, 102 

Lagerros, J. S. V. 1996, Astronomy and Astrophysics, 310, 1011 

Lagen'os, J. S. V. 1997, Astronomy and Astrophysics, 325, 1226 

Lagerros, J. S. V. 1998, Astronomy and Astrophysics, 332, 1123 

Larson, S., Beshore, E., Hill, R., et al. 2003, in Bulletin of the American Astronomical Society, Vol. 35, AAS/Division for Planetaiy Sciences 
Meeting Abstracts #35, 982 

Lauretta, D. S., Barucci, M. A., Bierhaus, E. B., et al. 2012, LPI Contributions, 1667, 6291 

Lupishko, D. R, Velichko, R R, Kazakov, V. V, & Shevchenko, V. G. 1987, Kinematika i Fizika Nebesnykh Tel, 3, 92 
Mainzer, A., Bauer, J., Grav, T., et al. 2011a, Astrophysical Journal, 731, 53 


30 


Mainzer, A., Bauer, J., Grav, T., et al. 2014, Astrophysical Journal, 784, 110 
Mainzer, A., Grav, T., Bauer, J., et al. 2011b, Astrophysical Journal, 743, 156 
Mainzer, A., Grav, T., Masiero, J., et al. 2011c, Astrophysical Journal, 736, 100 

Marchi, S., Delbo’, M., Morbidelli, A., Paolicchi, P., & Lazzaiin, M. 2009, Monthly Notices of the Royal Astronomical Society, 400, 147 
Marchis, F., Durech, J., Castillo-Rogez, J., et al. 2014, Astrophysical Journal, Letters, 783, L37 
Marchis, F., Kaasalainen, M., Horn, E. F. Y, et al. 2006, Icarus, 185, 39 

Marciniak, A., Michalowski, T., Hirsch, R., et al. 2009, Astronomy and Astrophysics, 508, 1503 
Masiero, J. R., Mainzer, A. K., Grav, T., et al. 2011, Astrophysical Journal, 741, 68 
Masiero, J. R., Mainzer, A. K., Grav, T., et al. 2012, Astrophysical Journal, Letters, 759, L8 
Mellon, M. T., Jakosky, B. M., Kielfer, H. H., & Christensen, P. R. 2000, Icarus, 148, 437 
Mueller, M. 2012, ArXiv e-prints 

Mueller, M., Marchis, F., Emery, J. P., et al. 2010, Icarus, 205, 505 
Muinonen, K., Belskaya, 1. N., Cellino, A., et al. 2010, Icarus, 209, 542 

Muller, T. G., Hasegawa, S., & Usui, F. 2014, Publications of the Astronomical Society of Japan, 66, 52 
Muller, T. G., O’Rourke, L., Barucci, A. M., et al. 2012, Astronomy and Astrophysics, 548, A36 

Okada, T., Fukuhara, T., Tanaka, S., et al. 2014, in Lunar and Planetary Inst. Technical Report, Vol. 45, Lunai' and Planetary Science Conference, 
1201 

Oszkiewicz, D. A., Muinonen, K., Bowell, E., et al. 2011, Journal of Quantitative Spectroscopy & Radiative Transfer, 112, 1919 
Piironen, J., Lagerkvist, C., Torppa, J., Kaasalainen, M., & Warner, B. 2001, in Bulletin of the American Astronomical Society, Vol. 33, Bulletin 
of the American Astronomical Society, 1562 
Pilcher, R, Benishek, V., Briggs, J. W., et al. 2012, Minor Planet Bulletin, 39, 141 
Pravec, P, Sai'ounova, L., & Wolf, M. 1996, Icarus, 124, 471 

Press, W. H., Elannery, B. P, & Teukolsky, S. A. 1986, Numerical recipes. The art of scientific computing (Cambridge: University Press, 1986) 
Rozitis, B. & Green, S. F. 2014, Astronomy and Astrophysics, 568, A43 
Rozitis, B., MacLennan, E., & Emery, J. P 2014, Nature, 512, 174 
Salisbury, J. W., D’Aria, D. M., & Jarosewich, E. 1991, Icarus, 92, 280 

S^neczky, K., Szabo, G., & Kiss, L. L. 1999, Astronomy and Astrophysics Supplement Series, 137, 363 
Skiff, B. A., Bowell, E., Koehn, B. W., et al. 2012, Minor Planet Bulletin, 39, 111 
Spencer, J. R. 1990, Icaims, 83, 27 

Spencer, J. R., Lebofsky, L. A., & Sykes, M. V. 1989, Icaims, 78, 337 
Stephens, R. D. 2009a, Minor Planet Bulletin, 36, 59 
Stephens, R. D. 2009b, Minor Planet Bulletin, 36, 18 

Szabo, G. M., Cs^, B., Sameczky, K., & Kiss, L. L. 2001, Astronomy and Astrophysics, 375, 285 
Tedesco, E. R, Noah, P. V., Noah, M., & Price, S. D. 2002, Astronomical Journal, 123, 1056 
Tholen, D. J. 1984, PhD thesis, Arizona Univ., Tucson. 

Tholen, D. J. 1989, in Asteroids II, ed. R. P. Binzel, T. Gehrels, & M. S. Matthews, 1139-1150 

Usui, R, Kuroda, D., Muller, T. G., et al. 2011, Publications of the Astronomical Society of Japan, 63, 1117 

Velichko, F. P, Krugly, Y. N., Lupishko, D. R, & Mokhamed, R. A. 1990, Astronomicheskij Tsirkulyar, 1546, 39 

Velichko, F. P, Psarev, V. A., Kiselev, N. N., et al. 2013, in Lunar and Planetary Science Conference, Vol. 44, Lunar and Planetary Science 
Conference, 2372 

Vernazza, P, Delbo, M., King, P. L., et al. 2012, Icarus, 221, 1162 

Warner, B. 2000, Minor Planet Bulletin, 27, 4 

Warner, B. D., Harris, A. W., & Pravec, P. 2009, Icarus, 202, 134 


31 


Wisniewski, W. Z., Michaiowski, T. M., Hams, A. W., & McMillan, R. S. 1997, Icarus, 126, 395 
Wright, E. L., Eisenhardt, R R. M., Mainzer, A. K., et al. 2010, Astronomical Journal, 140, 1868 
Yoshida, R, Dermawan, B., Ito, T., et al. 2004, Publications of the Astronomical Society of Japan, 56, 1105 
Zimbelman, J. R. 1986, Icarus, 68, 366 


32 


Appendix A. WISE data 


A.l. Our estimation of the uncertainties of WISE thermal infrared data 






Figure 8: Distribution of normalized differences f \2 for all pairs of close points in our WISE thermal data set in filters W3 (left panel) and W4, 
respectively (right panel), and the best normal-fit approximations (lines). 


The purpose of this section is to describe the method we used to independently estimate the accuracy of the WISE 
thermal infrared data. Our estimate is based on the examination of WHAT we call double detections of asteroids. Due 
to an overlap of about 10% between the areas of the field of view of two subsequently observed frames by WISE, 
we can sometimes find two flux measurements of an asteroid separated by the ~11 seconds, which is the cadence 
between the WISE frames. We have identified ~400 such double detections in the WISE data set for asteroids with 
convex models. The change in the thermal flux due to orbital and rotational evolution during this short time interval 
can be neglected, because it is smaller than the rotational spin-barrier period (~2.2 h) by a factor of ~700. The flux of 
an asteroid with a double peaked thermal lightcurve and 2.2 hour rotational period changes only by ~0.03% in 11 s, 
which is significantly smaller than the typical flux error of ~l-2%. Thus we can use these measurements to estimate 
the errors. 

Let us assume two measurements separated by 11 s with fluxes /i and /2 (/i is always the first one observed) and 
reported uncertainties cri and cr 2 . We also assume that these independent points are drawn from normal distributions 
N(//, cTj) and N(yu, cr^), respectively (the mean fi is unknown). The difference fi — f 2 is also a random value from a 
normal distribution N(0, cr^ 2 ^, where dispersion cri 2 corresponds to 


cr 12 



-HO-2. 


(3) 


By normalizing the difference /i - / 2 by cri 2 , we get a random value /n from a distribution N(0,1). We calculate 
/i 2 for all double detections we find in Alters W3 and W4 separately, and test the hypothesis that the fluxes /i and 
/2 are drawn from normal distributions N(yr, (Tj) and N(jU, cr^). If we obtain values /n following distribution N(0,1), 
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the uncertainties <ti and 0-2 correspond to standard Icr values. However, by fitting these distributions of normalized 
differences /n, we determine that they correspond to N(0.19 + 0.07, (1.40 + 0.07)‘) and N(0.02 + 0.06, (1.31 + 0.06)^), 
respectively (see Fig. [^. The mean for W4 data is close to zero, and thus the data are not significantly offset. On 
the other hand, there is a small offset in the W3 filter consistent with the flux measurements from one position on the 
frame (i.e., /i) to be systematically higher than the fluxes from the position on the opposite side of the frame (this 
is given by the particular scanning law). The reason could be vignetting and may not apply for the central parts of 
the frame. Our analysis shows that more realistic error bars in filters W3 and W4 are a factor ~1.4 and ~1.3 larger 
than the ones reported from the wise catalog. We also use the Kolmogorov-Smirnov test to test two zero hypotheses: 
(i) values /12 are drawn from a normal distribution N(0,1), and (ii) values /12 divided by the factors ~1.4 and ~1.3 are 
drawn from a normal distribution N(0,1). The hypothesis (i) is rejected at level 10“^ for W3 filter and 0.006 for W4 
filter. On the other hand, hypothesis (ii) cannot be rejected at level 0.92 and 0.87. 


A.2. Cross-calibration of the WISE data against Spitzer 



Figure 9: Chi-square recalculation for different offsets (±10%) introduced to W4 thermal data (W3 data are fixed) of asteroids (771) Libera and 
(2606) Odessa. The red lines correspond to the curves obtained for the original thennal data. It shows the sensitivity of the TPM fits to the possible 


offsets observed by 


Jarrett et al. 


201 ih in their WISE-Spitzer cross-calibration study. Note that these offsets were characterized for stars and may 


not be applicable to asteroids. 


As already introduced in Sect l2.2l 


.Tarrett et al 


(1201 ih studied the accuracy of the absolute calibration of WISE 


data by performing an analysis of WISE-Spitzer flux cross-calibration of a number of calibration stars and one galaxy 
situated near the poles of the ecliptic. They found that (i) the photometry for individual objects is stable for the 
whole cryogenic phase within less than 1%, but (ii) there is an rms scatter around the zero level of 4.5% and 5.7% of 
the WISE zero magnitudes in filters W3 and W4 when several objects are examined (the offset is different for each 
individual star). 

Given that filters W3 and W4 are independent, we expect the offsets in these filters to be different for each 
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individual asteroid. As a result, the fluxes will be changed by different multiplicative factors and the temperature will 
be different. To qualitatively estimate the corresponding change in the best-fitting values of parameters derived from 
the TPM, we perform the following test: we fix the observed fluxes in Alter W3 and introduce various multiplicative 
factors in W4 data spanning +10% with a step of 2% and run the TPM scheme with the fixed shape models. We 
note that +10% is is a slightly overestimated range, because the rms scatter of ^5% in both Alters should be summed 
in quadrature, which results in =s7%. In Fig. |9] we show the thermal inertia curves for different offsets introduced 
to W4 data of asteroids (771) Libera and (2606) Odessa. The effect on the thermal inertia is only apparent for the 
largest offsets, and, in general, the solutions remain well within the error bars given in Table |4] Indeed, for the most 
extreme cases, the best-fitting thermal inertia F could change by up to 50% (which is not outside the limits for error 
bars reported in literature), the size D by up to 10% and the geometric visible albedo pv by up to 20%. This range of 
uncertainties represents the extreme cases that would be introduced by the discrepancy seen in the cross-calibration 
of WISE and Spitzer. Note that changing the fluxes by the same multiplicative factor will not affect the value of the 
thermal inertia or the surface roughness, and that it would only affect the diameter. 


Appendix B. Rotational phase of the shape model 

In this section, we justify in more detail the necessity to optimize the rotational phase (po in the TPM. 

All shape models we use here are derived by the lightcurve inversion technique. This gradient-based method 
searches a large parameter space that includes sidereal rotational period, pole orientation, shape and scattering pa¬ 
rameters, and converges to all local minima and essentially finds the deepest (global) minimum. The difference AP 
between two local minima in the parameter space of rotational periods corresponds to 


^P ^ — , 
2T 


(4) 


where P is the sidereal rotational period and T = 7’2 - Ti is the timespan of the photometric data (IKaasalainen et al 


200 ih and T\ and T 2 are the epochs of the first and last observations of the asteroid within the photometric data set. 
The meaning of Eq. (HI is the following: the rotational phase shift during time T due to a change in period of AP is 
180° (assuming double-peaked sinusoidal lightcurve), thus the difference between the corresponding maxima of the 
lightcurves with P and P + AP after time T is exactly 180°. 

Let us assume that the uncertainty of the initial shape model orientation d(po due to the noise in the optical 
lightcurves corresponds to ~10 degrees, thus we can write for the uncertainty of the rotational period dP 



(5) 


This means that the orientation of the shape model is only known to this level of certainty within the coverage of the 
optical observations. Because the WISE data were acquired outside the timespan T (Tw > T 2 ), the uncertainty in the 
shape model’s initial orientation 6(po is even larger: 
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50o[deg] = 10 + 360 ^'^ 


(6) 


So, values of the uncertainty dtpo of the rotational phase are always at least 10° purely because of the uncertainty in 
P. Some of the shape models from the DAMIT database we use here were derived more than 10 years before WISE 
observations, thus the values of 6(po are large. The most extreme case is the asteroid (1036) Ganymed, for which 
d(p() ~ 540°, thus it is completely unconstrained (the model is based on the data from only one apparition in 1985). 
The expected uncertainties 6(po for other studied asteroids are < 20°. Of course, d(po is dependent on the value of 6P, 
for which we assume a reasonable value. 
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Figure 10: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for asteroid (624) Hektor (revised model). 
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Figure 11: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for asteroid (771) Libera (DAMIT model). There are two groups observations separated by few months (top and bottom 
panels). The offset between several model points and the synthetic lightcurve in the lower panels is caused by the fact that the obseiwations span four 
days. In this time, the thermal lighcurve changes because of the variation in the observing geometry. However, the synthetic lightcurve corresponds 
to the beginning of this interval. 
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Figure 12: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for both pole solutions of asteroid (832) Karin (DAMIT model). The small offset between several model points and the 
synthetic lightcurve is caused by the fact that the observations span several rotational periods. In this time, the thermal lighcurve changes because 
of the variation in the observing geometry. However, the synthetic lightcurve corresponds to the beginning of this interval. 
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Figure 13: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for asteroid (1036) Ganymed (revised model). There are two groups observations separated by few months (top and 
bottom panels). The small offset between several model points and the synthetic lightcurve is caused by the fact that the observations span several 
rotational periods. In this time, the thermal lighcurve changes because of the variation in the observing geometry. However, the synthetic lightcurve 
corresponds to the beginning of this interval. 
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Figure 14: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for both pole solutions of asteroid (1472) Muonio (DAMIT model). 
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Figure 15: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for asteroid (1627) Ivar (revised model). 
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Figure 16: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for asteroid (1865) Cerberus (revised model). The small offset between several model points and the synthetic lightcurve 
is caused by the fact that the observations span several rotational periods. In this time, the thermal lighcurve changes because of the variation in the 
observing geometry. However, the synthetic lightcurve corresponds to the beginning of this interval. 
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Figure 17: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for asteroid (1980) Tezcatlipoca (revised model). There are two groups observations separated by few months (top and 
bottom panels). 
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Figure 18: Comparison between the observed (WISE) and modeled (best TPM fit) thermal IR fluxes in filters W3 and W4 and synthetic thermal 
lightcurves (dotted lines) for both pole solutions of asteroid (2606) Odessa (DAMIT model). 
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